*------------------------------------------------------------------------------*
*------------------------------------------------------------------------------*
*------------------------------------------------------------------------------*

* Authors: Ranjit Lall & David Vilalta

* Paper: Varieties of Anxieties: Disaggregating Emotion and Voting Behavior in the COVID-19 Era

* What for: Replicates Tables and Figures For CIS Section

*------------------------------------------------------------------------------*
*------------------------------------------------------------------------------*
*------------------------------------------------------------------------------*

clear all 

* Insert folder path
global folder ".../Insert folder path"

cd "${folder}"

*Local folders
global do_files "Do_Files/"
global dta "Final_data/"
global tables "Results/Tables/"
global figures "Results/Figures/"


use ${dta}Finaldata_CIS, clear

global parties proconf anticonf pod psoe cs pp vox

global vars_3 covideffectsb

global parties2 proconf anticonf

********************************************************************************
* FIGURE 2 / FIGURE A3
********************************************************************************


foreach party in $parties{

eststo `party'_comb_1_hw: logit `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)

eststo `party'_comb_1_w: logit `party' worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)
		
eststo `party'_comb_1b: logit `party' covideffectsb worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)

eststo `party'_comb_1_ols_hw: reg `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)
		
eststo `party'_comb_1_ols_w: reg `party' worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)

eststo `party'_comb_1_olsb: reg `party' covideffectsb worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov)
	
}



* FIGURE 2 - TOP PANEL
					
coefplot  ///
		(pod_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pod_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(psoe_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(psoe_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(proconf_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(proconf_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(anticonf_comb_1_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(anticonf_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)),  ///
				keep(worr) eform  grid(none) ///
				vertical ////
				xlab(0.6 `" "(1)" "Podemos" "' 0.75 `" "(2)" "PSOE" "' 0.87 `" "(3)" "Ciudadanos" "' 0.99 `" "(4)" "PP" "' 1.13 `" "(5)" "Vox" "' 1.25 `" "(6)" "Pro-Lockdown" "' 1.38  `" "(7)" "Anti-Lockdown" "') ///
				yline(1, lcolor(gs10) lpattern(dash)) ///
				xline(0.67, lcolor(gs15) lpattern(solid)) ///
				xline(0.81, lcolor(gs15) lpattern(solid)) ///		
				xline(0.93, lcolor(gs15) lpattern(solid)) ///
				xline(1.06, lcolor(gs15) lpattern(solid)) ///
				xline(1.19, lcolor(gs10) lpattern(dash)) ///		
				xline(1.32, lcolor(gs15) lpattern(solid)) ///
				legend(order(1 "99% CI" 2 "95% CI" 3 "90% CI" 4 "Baseline" 8 "Controlling for Other Anxiety Measure") rows(1)) ///
				graphregion(color(white)) ///
				ytitle("Explanatory Variable:" "General COVID-19 Anxiety", margin(medium) box fcolor(gs13) lcolor(gs8) lwidth(medium) size(medium) width(100))  ///
				xtitle("Voting Intention (Political Party)") ///
				xsize(6) ysize(2)


* FIGURE 2 - BOTTOM PANEL
					
coefplot  ///
		(pod_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pod_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(psoe_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(psoe_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(proconf_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(proconf_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(anticonf_comb_1_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(anticonf_comb_1b, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)),  ///
				keep(covideffectsb) eform  grid(none) ///
				vertical ////
				xlab(0.6 `" "(1)" "Podemos" "' 0.75 `" "(2)" "PSOE" "' 0.87 `" "(3)" "Ciudadanos" "' 0.99 `" "(4)" "PP" "' 1.13 `" "(5)" "Vox" "' 1.25 `" "(6)" "Pro-Lockdown" "' 1.38  `" "(7)" "Anti-Lockdown" "') ///
				yline(1, lcolor(gs10) lpattern(dash)) ///
				xline(0.67, lcolor(gs15) lpattern(solid)) ///
				xline(0.81, lcolor(gs15) lpattern(solid)) ///		
				xline(0.93, lcolor(gs15) lpattern(solid)) ///
				xline(1.06, lcolor(gs15) lpattern(solid)) ///
				xline(1.19, lcolor(gs10) lpattern(dash)) ///		
				xline(1.32, lcolor(gs15) lpattern(solid)) ///
				legend(order(1 "99% CI" 2 "95% CI" 3 "90% CI" 4 "Baseline" 8 "Controlling for Other Anxiety Measure") rows(1)) ///
				graphregion(color(white)) ///
				ytitle("Explanatory Variable:" "Health-Weighted COVID-19 Anxiety", margin(medium) box fcolor(gs13) lcolor(gs8) lwidth(medium) size(medium) width(100))  ///
				xtitle("Voting Intention (Political Party)") ///
				xsize(6) ysize(2)


* FIGURE A3 - TOP PANEL

coefplot  ///
		(pod_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pod_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(psoe_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(psoe_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(proconf_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(proconf_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(anticonf_comb_1_ols_w, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(anticonf_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)),  ///
				keep(worr) eform  grid(none) ///
				vertical ////
				xlab(0.6 `" "(1)" "Podemos" "' 0.75 `" "(2)" "PSOE" "' 0.87 `" "(3)" "Ciudadanos" "' 0.99 `" "(4)" "PP" "' 1.13 `" "(5)" "Vox" "' 1.25 `" "(6)" "Pro-Lockdown" "' 1.38  `" "(7)" "Anti-Lockdown" "') ///
				yline(1, lcolor(gs10) lpattern(dash)) ///
				xline(0.67, lcolor(gs15) lpattern(solid)) ///
				xline(0.81, lcolor(gs15) lpattern(solid)) ///		
				xline(0.93, lcolor(gs15) lpattern(solid)) ///
				xline(1.06, lcolor(gs15) lpattern(solid)) ///
				xline(1.19, lcolor(gs10) lpattern(dash)) ///		
				xline(1.32, lcolor(gs15) lpattern(solid)) ///
				legend(order(1 "99% CI" 2 "95% CI" 3 "90% CI" 4 "Baseline" 8 "Controlling for Other Anxiety Measure") rows(1)) ///
				graphregion(color(white)) ///
				ytitle("Explanatory Variable:" "General COVID-19 Anxiety", margin(medium) box fcolor(gs13) lcolor(gs8) lwidth(medium) size(medium) width(100))  ///
				xtitle("Voting Intention (Political Party)") ///
				xsize(6) ysize(2)

* FIGURE A3 - BOTTOM PANEL
				
coefplot  ///
		(pod_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pod_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(psoe_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(psoe_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(proconf_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(proconf_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(anticonf_comb_1_ols_hw, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(anticonf_comb_1_olsb, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)),  ///
				keep(covideffectsb) eform  grid(none) ///
				vertical ////
				xlab(0.6 `" "(1)" "Podemos" "' 0.75 `" "(2)" "PSOE" "' 0.87 `" "(3)" "Ciudadanos" "' 0.99 `" "(4)" "PP" "' 1.13 `" "(5)" "Vox" "' 1.25 `" "(6)" "Pro-Lockdown" "' 1.38  `" "(7)" "Anti-Lockdown" "') ///
				yline(1, lcolor(gs10) lpattern(dash)) ///
				xline(0.67, lcolor(gs15) lpattern(solid)) ///
				xline(0.81, lcolor(gs15) lpattern(solid)) ///		
				xline(0.93, lcolor(gs15) lpattern(solid)) ///
				xline(1.06, lcolor(gs15) lpattern(solid)) ///
				xline(1.19, lcolor(gs10) lpattern(dash)) ///		
				xline(1.32, lcolor(gs15) lpattern(solid)) ///
				legend(order(1 "99% CI" 2 "95% CI" 3 "90% CI" 4 "Baseline" 8 "Controlling for Other Anxiety Measure") rows(1)) ///
				graphregion(color(white)) ///
				ytitle("Explanatory Variable:" "Health-Weighted COVID-19 Anxiety", margin(medium) box fcolor(gs13) lcolor(gs8) lwidth(medium) size(medium) width(100))  ///
				xtitle("Voting Intention (Political Party)") ///
				xsize(6) ysize(2)


********************************************************************************
* FIGURE A2
********************************************************************************
	   
histogram worr, percent discrete color(navy) lcolor(white)  width(0.5) ///
    xlabel(1 "None" 2 "Average" 3 "Not Much" 4 "Quite A Bit" 5 "A Lot") ///
    graphregion(fcolor(white)) ///
	xtitle("Response Options") ///
	ytitle("Share of Respondents (%)")
	   

********************************************************************************
* FIGURE 3
********************************************************************************

* Age

foreach n of numlist 1/6 {

eststo anx_age_`n'_1: reghdfe covideffectsb i.jobtype_old i`n'.agegroup sex i.class left_right_scale i.studies i.laboursit ln_num_casos_p_P, vce(cluster prov) a(i.prov i.date2)

}

	
* Class

foreach n of numlist 1/5 {

eststo anx_class_`n'_1: reghdfe covideffectsb i.jobtype_old i.agegroup sex i`n'.class left_right_scale i.studies i.laboursit ln_num_casos_p_P, vce(cluster prov) a(i.prov i.date2)

}


* Studies

foreach n of numlist 1/4 {

eststo anx_studies_`n'_1: reghdfe covideffectsb i.jobtype_old i.agegroup sex i.class left_right_scale i`n'.studies i.laboursit ln_num_casos_p_P, vce(cluster prov) a(i.prov i.date2)

}


* Job type

foreach n of numlist 1/11 {

eststo anx_jobtype_`n'_1: reghdfe covideffectsb i`n'.jobtype_old i.agegroup sex i.class left_right_scale i.studies i.laboursit ln_num_casos_p_P, vce(cluster prov) a(i.prov i.date2)

}


* Labour situation
foreach n of numlist 1/4 {

eststo anx_laboursit_`n'_1: reghdfe covideffectsb i.jobtype_old i.agegroup sex i.class left_right_scale i.studies i`n'.laboursit ln_num_casos_p_P, vce(cluster prov) a(i.prov i.date2)

}


coefplot ///
		(anx_age_1_1, keep(1.agegroup))  ///
		(anx_age_2_1, keep(2.agegroup))  ////
		(anx_age_3_1, keep(3.agegroup))  ///
		(anx_age_4_1, keep(4.agegroup))  ///
		(anx_age_5_1, keep(5.agegroup))  ///
		(anx_age_6_1, keep(6.agegroup))  ///
		(anx_class_1_1, keep(1.class))  ///
		(anx_class_2_1, keep(2.class))  ///
		(anx_class_3_1, keep(3.class))  ///
		(anx_class_4_1, keep(4.class))  ///
		(anx_class_5_1, keep(5.class))  ///		
		(anx_studies_1_1, keep(1.studies))  ///
		(anx_studies_2_1, keep(2.studies))  ///
		(anx_studies_3_1, keep(3.studies))  ///
		(anx_studies_4_1, keep(4.studies))  ///
		(anx_laboursit_1_1, keep(1.laboursit))  ///
		(anx_laboursit_2_1, keep(2.laboursit))  ///
		(anx_laboursit_3_1, keep(3.laboursit))  ///
		(anx_laboursit_4_1, keep(4.laboursit)) ///
		(anx_jobtype_1_1, keep(1.jobtype_old))  ///
		(anx_jobtype_2_1, keep(2.jobtype_old))  ///
		(anx_jobtype_5_1, keep(5.jobtype_old))  ///
		(anx_jobtype_6_1, keep(6.jobtype_old)), ///
		drop(5.laboursit 6.laboursit ln_num_casos_p_P _cons sex) ///
		coeflabels (1.laboursit = "Working" 2.laboursit = "Retired" 3.laboursit = "Unemployed" 4.laboursit = "Student" ///
		1.studies = "No Education" 2.studies = "Primary Education" 3.studies = "         Secondary Education" 4.studies = "Higher Education" ///
		1.class = "Working Class" 2.class = "  Lower Middle Class" 3.class = "Middle Class" 4.class = "  Upper Middle Class" 5.class = "Upper Class" ///
		1.agegroup = "18-24" 2.agegroup = "25-34" 3.agegroup = "35-44" 4.agegroup = "45-54" 5.agegroup = "55-64" 6.agegroup = "65+" ///
		1.jobtype_old = "Manager" 2.jobtype_old = "Scientist/Intellectual" 3.jobtype_old = "Mid-Level Professional" 4.jobtype_old = "Administrator" 5.jobtype_old = "Service Worker" 6.jobtype_old = "Agricultural Worker" 7.jobtype_old = "Office Worker" 8.jobtype_old = "Machine Operator" 9.jobtype_old = "Elementary Worker" 10.jobtype_old = "Security Worker")  ///
		 graphregion(color(white)) ///
		headings(1.laboursit = "{bf:Labor Situation}"  ///
              1.studies = "{bf:  Education Level}" ///
			  1.class = "{bf:Social Class}" ///
			  1.agegroup = "{bf:Age Group}" ///
			  1.jobtype_old = "{bf:Job Type}") ///                       
		xline(0, lcolor(gs10) lpattern(dash))  norecycle  ///
		legend(order(1 "99% CI" 2 "95% CI" 3 "90% CI") rows(1) region(lstyle(none)) symxsize(4)) ///
		xtitle(Coefficient Estimate) ///  
		subtitle(, fcolor(gs13) lcolor(gs8) margin(medium) lwidth(medium) ) ///
		ysize(13) xsize(8) graphregion(margin(large)) ///
		mcolor(navy) msize(medium) msymbol(t) ciopts(color(navy teal ltblue)) level(99 95 90)
		
********************************************************************************
* FIGURE A4
********************************************************************************


foreach party in $parties{

eststo `party'_worr_1r: logit `party' worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.date, vce(cluster prov), if covideffectsb!=.
		

eststo `party'_worr_2r: logit `party' worr i.jobtype_old i.agegroup sex i.class i.studies i.laboursit `party'prev i.prov ln_num_casos_p_P i.prov i.ccaa##i.date, vce(cluster prov), if covideffectsb!=.

}

					
coefplot  ///
		(pod_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pod_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(psoe_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(psoe_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(cs_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(pp_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(vox_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(proconf_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(proconf_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)) ///
		(anticonf_worr_1r, mcolor(ebblue) msize(large) msymbol(t) ciopts(color(ebblue midgreen cranberry)) level(99 95 90))  ///
		(anticonf_worr_2r, mcolor(ebblue) msize(large) msymbol(s) ciopts(color(ebblue midgreen cranberry)) level(99 95 90)),  ///
				keep(worr) eform  grid(none) ///
				vertical ////
				xlab(0.6 `" "(1)" "Podemos" "' 0.75 `" "(2)" "PSOE" "' 0.87 `" "(3)" "Ciudadanos" "' 0.99 `" "(4)" "PP" "' 1.13 `" "(5)" "Vox" "' 1.25 `" "(6)" "Pro-Lockdown" "' 1.38  `" "(7)" "Anti-Lockdown" "') ///
				yline(1, lcolor(gs10) lpattern(dash)) ///
				xline(0.67, lcolor(gs15) lpattern(solid)) ///
				xline(0.81, lcolor(gs15) lpattern(solid)) ///		
				xline(0.93, lcolor(gs15) lpattern(solid)) ///
				xline(1.06, lcolor(gs15) lpattern(solid)) ///
				xline(1.19, lcolor(gs10) lpattern(dash)) ///		
				xline(1.32, lcolor(gs15) lpattern(solid)) ///
				legend(off) ///
				graphregion(color(white)) ///
				ytitle("Treatment: General COVID-19 "" Anxiety" " ", margin(medium) box fcolor(gs13) lcolor(gs8) lwidth(medium) size(medium) width(100))  ///
				xsize(6) ysize(2)

********************************************************************************		
* TABLE A5
********************************************************************************

foreach party in $parties2{

logit `party' covideffectsb  i.prov i.date, vce(cluster prov)

logit `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit i.prov i.date, vce(cluster prov)
		
logit `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit ln_num_casos_p_P  i.prov i.date, vce(cluster prov)

logit `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit ln_num_casos_p_P i.party_sympathy_num i.prov i.date, vce(cluster prov)

logit `party' covideffectsb i.jobtype_old i.agegroup sex i.class i.studies i.laboursit ln_num_casos_p_P left_right_scale i.prov i.date, vce(cluster prov)

}

		
********************************************************************************
* TABLE A11
********************************************************************************		
				  
use ${dta}Finaldata_CIS, clear
				  
* PANEL A

reghdfe worr ln_num_casos_p_P, vce(cluster prov) a(i.prov i.ccaa##i.date2)

reghdfe worr ln_num_casos_p_P i.agegroup sex i.class left_right_scale i.studies, vce(cluster prov) a(i.prov i.ccaa##i.date2)

reghdfe worr ln_num_casos_p_P, vce(cluster ccaa) a(i.prov i.ccaa##i.date2)

reghdfe worr ln_num_casos_p_P i.agegroup sex i.class left_right_scale i.studies, vce(cluster ccaa) a(i.prov i.ccaa##i.date2)

reghdfe worr ln_num_casos_p_P, vce(cluster ccaa_date) a(i.prov i.ccaa##i.date2)

reghdfe worr ln_num_casos_p_P i.agegroup sex i.class left_right_scale i.studies, vce(cluster ccaa_date) a(i.prov i.ccaa##i.date2)

		
* PANEL B	

reghdfe worr symptoms, vce(cluster prov) a(i.prov i.ccaa##i.date2)

reghdfe worr symptoms i.agegroup sex left_right_scale i.class i.studies, vce(cluster prov) a(i.prov i.ccaa##i.date2)

reghdfe worr symptoms, vce(cluster ccaa) a(i.prov i.ccaa##i.date2)

reghdfe worr symptoms i.agegroup sex left_right_scale i.class i.studies, vce(cluster ccaa) a(i.prov i.ccaa##i.date2)

reghdfe worr symptoms, vce(cluster i.ccaa##i.date2) a(i.prov i.ccaa##i.date2)

reghdfe worr symptoms i.agegroup sex left_right_scale i.class i.studies, vce(cluster i.ccaa##i.date2) a(i.prov i.ccaa##i.date2)
	
		
		
		
		
		